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Abstract. We present a pilot study of atmospheres of accreting magnetic white dwarfs irradiated by intense fluxes at ultraviolet 
to infrared wavelengths. The model uses a standard LTE stellar atmosphere code which is expanded by introducing an angle- 
dependent external radiation source. The present results are obtained for an external source with the spectral shape of a 10000 K 
blackbody and a freely adjustable spectral flux. The model provides an explanation for the observed largely filled-up Lyman 
lines in the prototype polar AM Herculis during its high states. It also confirms the hypotheses (i) that irradiation by cyclotron 
radiation and other radiation sources is the principle cause for the large heated polar caps surrounding the accretion spots on 
white dwarfs in polars and (ii) that much of the reprocessed light appears in the far ultraviolet and not in the soft X-ray regime 
as suggested in the original simple theories. We also briefly discuss the role played by hard X-rays in heating the polar cap. 
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1. Introduction 


The white dwarfs in magnetic cataclysmic variables (mCVs, 
more specifically polars) accrete in restricted regions near their 
magnetic poles. In most of them, the infalling matter is heated 
in a free-standing shock above the surface of the white dwarf 
and cools by the emission of X-rays and cyclotron radiation. 
Early models assumed that the stand-off distance of the shock 
is small and much of the emission by the post-shock plasma is 
absorbed by the atmosphere below and immediately around the 
accretion spot. The implied irradiation of this small section of 
the atmosphere by intense fluxes of hard X-ray bremsstrahlung 
and cylotron radiation leads to an expected temperature of the 
heated atmosphere of kT =* 25 eV (T =* 3 x 10 5 K), to re¬ 
emission in the soft X-ray regime with a quasi-blackbody spec¬ 
trum, and to an energy balance = f\ L rvr + (1 - ,4 hx ) L hx ] 
Lamb & Mastersil979l). where the three 


trum, and to an energ 

dKing&Lasotall979t 


luminosity components refer to soft X-rays (sx), hard X-rays 
(hx), and cyclotron radiation (eye), the geometry factor / $ 
1 /2 accounts for the fraction of the luminosity intercepted by 
the white dwarf, and the ha rd X-ray albedo Ah, redu ces the 
efficiency of X-ray heating (Ivan Teeseling et al.l 1994] ). Later 
on, it was realized that the accretion region is hi ghly struc¬ 
tured dKuiiners & Pringlelll982t iFrank et afll 19881) with indi¬ 
vidual sections receiving vastly different mass flow densities 
and displaying substantially different emission characteristics. 
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The plasma in tenuous tall shocks cools preferentially by cy¬ 
clotron radiation, the plasma in dense low-lying shocks by hard 
X-ray bremsstrahlung, and for the highest mass flow densi¬ 
ties, the shocks are buried in the atmosphere and the primary 
bremsstrahlung is reprocessed into soft X-rays, a component 
physically different from but observationally difficult to dis¬ 
tinguish from the originally suggested soft X-ray blackbody 
component. This “soft X-ray puzzle” was the subject of an ex- 
tended debate over more than two decades (e.g., iBeuermamil 
l2004tfRamsav & Croppeil2004 and references therein). 


Quite unexpectedly, observational support for the repro¬ 
cessing scenario came from studies of the far ultraviolet 
spectra of the prototype polar AM Herculis, which showed 


the FUV (Heise & Verbunt 

1988; Gansicke et al. 1993 1998; 

Mauche & Raymond 1998[ 

Greeley et al. 1999; Gansicke et al. 

2005). The pronounced orbital modulation of the FUV flux 


indicates a large heated polar cap, which covers ~ 10% of 
the white dwarf surface. The cap reaches a peak temperature 
of ~ 3.5 x 10 4 K in low states when accretion nearly ceases 
and becomes much hotter in high states. Iciiinsicke et~al. ( 1995 ) 
demonstrated that the excess FUV flux of the polar cap quanti¬ 
tatively agrees with the sum of cyclotron and X-ray fluxes from 
AM Her, both in the high and the low states, thus confirming 
the energy balance predicted by the simple models, with the 
decisive difference, however, that the heated polar cap covers 
a much larger area than the proper accretion spot and the re¬ 
processed flux emerges in the FUV and not at soft X-ray wave- 
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lengths. AM Her is the best observed polar, but FUV studies of 
other polars suggest that the situation encountered in AM Her 
is quite typical of the class (e.g.. lArauio-Betancor et alJ2005l) . 

The dominant FUV spectral feature of white dwarfs in po¬ 
lars are the Lyman absorption lines, which are deep and broad 
in low states, mimicking the pure hydrogen spectra of DA white 


dwarfs, but become filled u 

p almost entirely in states of high 

accretion (Heise & Verbun 

1988; 

Gansicke et al. 1995, 19981 

Mauche & Raymond 1998E 

Greelev et al. 1999; Gansicke et al. 

2005). Gansicke et al. 

(1998) first tried to model the filled- 


up Lyman lines using a simple toy model and concluded that 
the temperature stratification in the white dwarf atmosphere is 
severely affected over the entire irradiated polar cap. 

A spot due to reprocessing covering 10% of the white 
dwarf surface requires a point-like source at a height of 0.25 
white dwarf radii or an extended source at a lower height. 
The best candidates are optical/infrared cyclotron and to some 
extent hard X-ray emission from tall shocks, supplemented 
by the FUV to optical emission from the pre-shock accre¬ 
tion stream. A further potential source is shockless accre¬ 
tion at very low mass flow densities, m < 10 3 g cm 2 s~'. 
The infalling particles do not directly heat the lower at¬ 
mosphere but create a hot corona and the atmosphere be¬ 
low is ag ain heated radiatively iThomson & CawthorneIT987l 
Woelk & Beuermanni 11992}) . The task is, therefore, to con¬ 


sider the properties of atmospheres heated by a large exter¬ 
nal flux of ultraviolet/optical/infrared radiation. Modeling the 
UV-emitting spots in polars requires us to consider incoming 
radiation fluxes that exceed the flux emerging from the un¬ 
heated atmosphere of the white dwarf by factors up to ~ 200 
in the high state. Irradiated stellar atmospheres have been stud¬ 
ied by several authors for a variety of conditions which differ, 
however, from those cons i dered here (e.g. , L 2 ndon_£DdJj - 2 £jl 
van Teeseling et all 1994 iBrett & Smithll 1 993t Barman et al. 
2004]). Tn this paper, we present a pilot study of the structure 
and the angle-dependent spectra of white dwarfs irradiated by 
an intense source of infrared to ultraviolet radiation. 

The paper is arranged as follows. In Sect. 2 we describe the 
stellar atmosphere code used to calculate the irradiated atmo¬ 
sphere models. Sect. 3 introduces the essential parameters of 
the accretion regions in polars, in Sect. 4 we report the results 
for simple irradiation geometries, which are applied to the pro¬ 
totype polar AM Her in Sect. 5. In Sect. 6, we discuss the limi¬ 
tations of the approach and lines for future research. 


2. Radiation transfer in irradiated atmospheres 

We assume a pure hydrogen elemental composition, as appro¬ 
priate for DA white dwarfs and for the sections of the atmo¬ 
sphere of accreting white dwarfs in mCVs that far away from 
the main accretion spot. We introduce the external flux as a 
boundary condition and neglect any other disturbance of the 
atmosphere as, e.g., an external pressure exerted by low-level 
accretion outside the main accretion spot. 

We use a standard LTE atmosphere code to solve the ra¬ 
diative transfer in a plane parallel geometry with an incoming 
flux /j n (ju. A) at optical depth r = 0, prescribed as a function 


of wavelength A and zenith angle # via the direction cosine 
/i = cos#. The radiative transfer equation 

dl( n r , ,, 

F ---= I(r,fi,A) - S(t,A) (1) 

OT 

accounts for absorption and isotropic electron scattering via the 
source function 

k cr 

S = - B +- J , (2) 

k + cr k + cr 

where k is the absorption coefficient, cr is the Thomson scatter¬ 
ing cross section, B the Planck function, J the angle-averaged 
intensity, and the dependencies on r and A have been sup¬ 
pressed. Line opacity is included in the form of the Lyman, 
Balmer, Paschen, and Brackett series of hydrogen, including 
excited levels up to « = 5. For the hot atmospheres consid¬ 
ered here, molecular absorption is unimportant. Although we 
apply the results to magnetic white dwarfs in CVs we neglect 
Zeeman splitting of the Lyman lines. This is appropriate for 
the moderate field strengths considered here (|Zf| SS 30 MG), 
because the Lyman lines are split by the normal Zeeman effect 
with Ad = eA 2 \B\/(4nm e c 2 ), where e is the elementary charge, 
m e the electron mass, and c the speed of light. Specifically, for 
the polar AM Her considered below, the observed field strength 
\B\ = 14 MG implies Ad = 10A for Lya, which is small com¬ 
pared with the 90A contributed by the Stark effect to the width 
of Ly a at the white dwarf effective temparature of 20000 K. 

The radiative transfer equation is solved by a standard 
iRvbickii d 1971 1> elimination scheme. To this end, the quantities 


w(r,/r, A) = [Kt,/i, A) + /(t, -fi,A)] (3) 

v(t,/i. A) = [I(t,/i, A) - /(t, -fi,A)] (4) 


are introduced for 0 < /./ < 1. With their help, Eq. Q is trans¬ 
formed to a second order differential equation 


2 d 2 u(T,/i,A) 
** dr 2 


u(t,/i , A) — S (r, /l) 


(5) 


which is solved numerically. At large optical depths, the diffu¬ 
sion approximation provides the standard boundary condition. 
The angle-dependent incident intensity h a (ji, A) is implemented 
by modifying the upper boundary condition 


1(0,-fi. A) = I m (ji, A) A 0 . 


( 6 ) 


For isotropic irradiation, 7j n is independent of /./, for irradiation 
at a single angle §0 with /jq -cos-&q, I m is described by a delta- 
function 6(ji - fin). The incoming and outgoing fluxes at r = 0 
are denoted by 

F in(0, A) = 2 7T f [/(0,-^,/i)] II dll . (7) 

Jo 

F out (0,A) = 2 7T f [1(0,/u. A)] \xd[x. (8) 

Jo 

Independent of the magnitude of the incoming flux, the integral 
quantities 

s-»oo /"»00 

F[ n = I F m (A)dA and F 0 m = FUA)dA (9) 

Jo Jo 
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Fig. 1. Spectral energy distribution at an angle of 6 — 80° rela¬ 
tive to the field direction for the post-shock emission region on 
a 0.6 Mr., white dwarf at a distance D = 10 pc with field strength 
\B\ = 30 MG, area A = 4x 10 16 cm 2 , and mass flow densities of 
m — 1 g cm^V 1 (dashed curve) and 0.1 g cm^V 1 (solid curve) 
(Fischer & Beuermannil20oH their Fig. 9). For comparison, a 
blackbody with temperature T = 10 4 K and the same emitting 
area is shown (dot-dashed curve). 


are related by F out = F m + F w a, where 


F wd — 



FMdA = crT* d 


( 10 ) 


is the flux emerging from the unheated white dwarf atmosphere 
with effective temperature T wd . 

The condition of radiative equilibrium. 



k(J -B)dA = 



8H 

X ~dr dA = ° ’ 


( 11 ) 


with FI = (1 /47r) F the Eddington flux, is achieved by appro¬ 
priate depth-dependent adjustment of the temperature. We use 
a correction method that is based on solving the momentum 
equations 


d\f K J ) = 

8t 2 


d(f K J) 

8t 


faJ 


( 12 ) 


with Eq. {U} as a constraint. Here, fy — H/J and fr = K/J 
are the Eddington factors with K = \ J_ I fj 2 d/i the scalar ra¬ 
diation pressure. A correction 6T to 74r) is obtained by substi¬ 
tuting B with 


B = B[T(t ) + 6T(t)] » B(T) + —^bT(r) (13) 

dT 

and solving Eq. (11 21 simultaneously with Eq. ill III for 5T. 


plasma, T e its electron temperature, and \B\ the field strength. 
The plasma temperature is of the order of 10 8 K at the shock 
and the plasma cools at nearly constant pressure ( n e oc 1 /7'J 
as it moves down. Since |2?| is also approximately constant 
throughout the post-shock region, the emissivities vary approx¬ 
imately as £brems TjT 3 ^ 2 and e cyc oc T®, i.e., bremsstrahlung 
emission is concentrated near the white dwarf surface, while 
cyclotron emission originates throughout the column and, on 
the average, high above the surface. For a shock-heated plasma 
dominated by cyclotron cooling, the stand-off distance of the 
shock is given by 

/j sh = 3.3 x 10 8 \B\- 2 - 6 M\- 2 cm, (14) 

where \B\-j is the field strength in units of 10 7 G and M\ is the 
white dwarf mass in solar units llFischer & BeuermannlbOoH 
their Eq. 21). Equation ill 4> remains valid as long as h s h stays 
small compared with the radius of the white dwarf. For a 
0.6 M 0 white dwarf with \B\ = 100MG, /z s |, amounts to a 
minute 4 x 10 5 cm, while for \B\ = 10 MG, Eq. (tl4l yields 
/z s h - 1.8 x 10 s cm or about 20% of the white dwarf radius. 
Large irradiated polar caps are, therefore, expected to exist 
preferentially in low-field mCVs. 

The angle-dependent spectra emitted by the accretion re¬ 
gions on mag netic white dwarfs have been calculated by 
iFischer & Beuermannl ( 2001 . see also Woelk & Beuermann, 
1996). Figure [T] shows an example of the emitted spectral en¬ 
ergy distribution in a direction almost perpendicular to the 
magnetic field. Cyclotron emission is largely optically thick 
in the infrared and assumes a quasi Rayleigh-Jeans spectrum 
with harmonic structure superimposed. The cyclotron spec¬ 
trum breaks off in the optical/UV, where the emission re¬ 
gion becomes optically thin. At high frequencies, the X-ray 
bremsstrahlung tail is seen, which contains negligible flux if 
the mass flow density m is sufficiently low. This is a conse¬ 
quence of the dependence of sterns on n 2 oc nr. For compari¬ 
son, we show a blackbody for a temperature of 10 4 K with the 
same emitting area. Naturally, its spectral flux in the infrared 
is lower by a factor of about 10 4 compared with the 10 8 K cy¬ 
clotron source. The comparison shows that the general spectral 
shape of the cyclotron emission mimics that of a blackbody of 
about 10 4 K. For the purpose of the present pilot study, the ex¬ 
act shape and angle dependence of the cyclotron emission is not 
critical. We, therefore, approximate the shape of the irradiation 
spectrum by a 10 4 K blackbody and scale the flux to match that 
of the cyclotron spectrum at the surface element in question. 
The reader should realize that the irradiation spectrum belongs 
to the plasma temperature of ~ 10 8 K and that the irradiation 
flux can, therefore, substantially exceed that of a 10 4 K black¬ 
body. 


4. Results 


3. Irradiated polar caps in AM Her stars 

The frequency-integrated bremsstrahlung and cyclotron emis¬ 
sivities vary as eb rem s oc n 2 T^ 2 and e cyc oc n e T e \B\ 2 , respec¬ 
tively, with » e the electron (and ion) number density of the 


All results presented in this paper have been obtained for irra¬ 
diation by a source that emits a spectrum with the shape of a 
blackbody of 7bb = 10 4 K. Test calculations for various val¬ 
ues of 7’bb between 3000 K and 25000 K suggest that the spec¬ 
trum emerging from the heated atmosphere is insensitive to 7’bb 
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A (A) 

Fig. 2. Example for the relation between intrinsic spectral flux 
/- w( ] of the unheated atmosphere of the white dwarf (dotted 
curve), the low-temperature irradiation F m (dashed curve), and 
the outgoing flux F out of the heated atmosphere (solid curve). 


as long as the fraction of the emitted energy shortward of the 
Lyman edge is small fa (ini liii). 

In a first step, we consider single surface elements irradi¬ 
ated (1) isotropically and (2) at a single zenith angle t). With re¬ 
spect to case 2, we note that our atmosphere code treats the in¬ 
tensity 7(r, 0) as azimuthally symmetric and, hence, interpretes 
the incoming intensity at r = 0 as given on a cone with open¬ 
ing angle i) rather than at a single direction //, <p with azimuth 
angle <p. As a consequence, the code neglects the azimuthal 
asymmetry in the outgoing intensity that occurs in nature for 
a single surface element irradiated unidirectionally from 0, i p. 
In a second step, we use case-2 models to calculate the inte¬ 
grated intensity from all surface elements of the white dwarf 
visible from a given direction 8, (f>, which indicates the direction 
towards the observer. We expect that the asymmetries present 
in the emissions from the individual elements largely cancel 
in the integral over all surface elements. With this caveat, we 
calculate the integrated spectrum received from an entire white 
dwarf, which possesses a large hot spot produced by either a 
point source or an extended source of total luminosity L lo¬ 
cated at a height h above the photosphere. 

We describe the wavelength-integrated irradiation flux F m 
of Eq. (9) in units of F wd of Eq. (10) by 

F in = xF wd (15) 

and consider irradiation parameters x up to 200. We use an 
unheated white dwarf atmosphere with effective temperature 
r wd = 2 x 10 4 K and logg = 8, the approximate paramet ers of 
the white dwarf in AM Her Jfiansicke et al.1179951 i2005 ). The 
chosen effective temperature i s typic al also of the white dwarfs 
in other polars llArauio-Betancor et al.ll2005l) . Note that with 
Ebb = 10 4 K and 7 wd = 2 x 10 4 K, F ln for x = 1 and x = 100 
amounts to 16 times and 1600 times the nominal blackbody 
flux, respectively. 


4.1. Case 1: Isotropic irradiation of a plane-parallel 
atmosphere 

In Figure [2] we show an example of the relation between the 
incoming and outgoing flux for x = 4.6, i.e., the outgoing flux 
F out = (1 +x) F wd = 5.6 F wd and the effective temperature of the 
heated surface element is 7’ eff = (1 + x) l/4 T wd = 3.08 x 10 4 K. 
The optical/infrared irradiation is reprocessed into the outgoing 
ultraviolet flux sufficiently high up in the atmosphere to cause 
an inversion in the temperature stratification and the formation 
of emission line cores. 

In Figure^ we show the emitted spectral flux and the tem¬ 
perature stratification for five selected values of F m = x F wd 
with x = 1.8,3.1,6.7,17.2 and 30.5. Our calculations cover a 
range of x — 0.1 to ~ 200, with a model for a lower value of 
x acting as a start model for the next higher x. For large x, in¬ 
stabilities in the temperature stratification develop (see e.g. the 
model for x = 30.5 in Fig. 0 right panel), which ultimately 
limit the manageable irradiation level. Since we are still far 
from the Eddington limit, we think that this deficiency can be 
mended in future applications. 

For moderate irradiation with x A 3, the temperature inver¬ 
sion (Fig.0 right panel) is barely noticeable and the tempera¬ 
ture stays essentially constant as a function of the Rosseland 
optical depth r ross until it merges into the temperature pro¬ 
file of the unheated white dwarf at large optical depths. For 
x > 3, a pronounced temperature inversion evolves and for still 
higher values, x <; 10, a hot corona develops at optical depths 
Tross ~ 10~ 2 . As a consequence, the outgoing spectra display 
filled-up Lyman lines for x $ 3 and emission cores beyond. 
The toy model devised by Ifiansicke et alJ 1 1998 1 for AM Her 
assumed a temperature independent of r ross in the outer atmo¬ 
sphere and yielded filled up lines very similar to the situation 
obtained here for x = 2 - 3. 

4.2. Case 2: Irradiation of a plane-parallel atmosphere 
at a fixed angle 

We have calculated a 2-D grid of model atmospheres for the 
same range of incident fluxes as above, x = 0.1 to 200, and 
for direction cosines yu = cos# = 0.1 to 1.0 in steps of 0.1, 
with // the angle under which the irradiation occurs. For each 
atmosphere with irradiation parameters x,/r, our database con¬ 
tains the intensities I JO, x,p) as a function of wavelength for 
a range of emission angles 8 , with cos 8 = 0.1 to 1.0 in steps 
of 0.1. Selected results for x = 1.0 and 6.7 and /i — 0.1 and 
1.0 are shown in Fig. 0 The left panel is for 8 = 0° (look¬ 
ing perpendicular into the irradiated atmosphere), the right one 
for 8 = 84° (looking under a grazing angle into the atmo¬ 
sphere). In both panels, the solid curves are for vertical irradia¬ 
tion 07 = 0°,/r = 1.0) and the dashed curves for grazing irradi¬ 
ation (# = 84°, n = 0.1). Finally, in both panels, the lower and 
upper pair of curves are for irradiation parameters x = 1.0 and 
6.7, respectively. The individual curves show that irradiation 
incident vertically on the atmosphere leads to weaker emission 
cores than grazing irradiation, which is reprocessed high up 
in the atmosphere. Radiation leaving the atmosphere vertically 
originates at larger depths and displays broader absorption lines 
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Fig. 3. Emerging spectral flux (left) and temperature stratification (right) for atmospheres isotropically irradiated with five differ¬ 
ent levels of the 10 4 K blackbody-like spectral flux with 1.8, 3.1, 6.7, 17.2, and 30.5 times the integrated flux /- w( i of the 2 x 10 4 K 
white dwarf atmosphere (solid curves, from bottom to top). Also shown is the spectral flux (left) and the temperature stratification 
(right) of the un-irradiated white dwarf (dotted curves). 




x (A) x (A) 

Fig. 4. Directional fluxes at emission angles 0 = 0° (looking perpendicular into the irradiated atmosphere, left panel) and 8 = 84° 
(looking under a grazing angle into the atmosphere, right panel) relative to the radial direction for irradiation at /.i = 1.0 (0 = 0°, 
solid curves) and /i = 0.1 ($ = 84°, dashed curves), shown for two levels of the incident flux, 1.0F w( j and 6.7 F m j (lower and 
upper set of curves). 


compared with the radiation leaving the atmosphere at large 8. 
Low-level irradiation (x = 1) preserves the limb darkening of 
the un-irradiated atmosphere, whereas grazing irradiation at the 
higher irradiation level (x = 6.7) leads to limb brightening in 
the continuum and even more so in the lines. While this is what 
one might expect intuitively, the quantitative results allow us 
to construct models of accretion-induced irradiation for more 
complicated geometries. 

4.3. External point source and extended source 

We consider two models of the radiation received from the en¬ 
tire white dwarf irradiated by a source of luminosity L. The 
simpler model assumes a point source located at a height h 
above the surface and the more involved one assumes an ex¬ 


tended source at the same height, consisting of an array of N 
point sources of equal luminosity L\ = L/N , with f = 1 to N. 
Fig-EJillustrates the irradiation geometry. 

For irradiation by a point source, the spectral flux ft ob¬ 
served at Earth is given by the sum over the contributions from 
the m surface elements visible for a given orientation of the 
spotted white dwarf from the direction of the observer, 

h = Jy2 Yj A J C0S IAd j' X J ’ M (16) 

7=1 

Here, D is the distance between the observer and the source, 
Aj and cos 8j the area and the cosine of the viewing angle of 
the y'-th surface element with 8j the angle between the normal 
to this element and the direction to the observer, and Xj and pj 
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Fig. 5. Left: Array of surface elements used in calculating models of white dwarfs irradiated by an extended source at a height 
h = 0.1 R v .d above the photosphere, shown for phases 0 = 0 (left) and 0 = 0.5 (right). Center: Definition of angles a and 0. 
Right: Irradiation of a single surface element by an extended source of area A at height h above the photosphere. 6 indicates the 
direction to the observer. 



Fig. 6. Effective temperature in the large irradiation-heated 
spot on a 0.6 M e white dwarf as a function of the angular dis¬ 
tance a of the irradiated surface element from the center of 
the source. Irradiation by a point source (solid curves) or an 
extended source of 0.3/? wd radius with J? wd = 8.6 x 10 8 cm 
(dashed curves) was considered, both located at a height of 
0.1 R w d- The luminosity is 10 32 for the lower curves and 
10 33 erg s-' for the upper curves. The horizontal dotted curve 
is for the unheated white dwarf. 


the irradiation parameters of the element as defined above, of 
which Xj is given by 


F inj = J_ L_ 

F w d F wd 47ni 2 


(17) 



Fig. 7. HST/GHRS spectrum of AM Her in the January 1997 
high state at orbital phases 0 = 0.1 (heated polar cap in front) 
and 0 = 0.5 (heated polar cap at horizon). The model spectra 
are for irradiation by an isotropically emitting source of 7.7 x 
10 32 erg s 1 with radius 0.3 R wd located at a height 0.1 /? wd and 
for a distance of AM Her of D - 79 pc. The contribution by 
the accretion stream (dashed curve) and the spectral flux of the 
unheated white dwarf (dotted curve, added to the dashed curve) 
are included for comparison. For the definition of the orbital 
phase see lfiansicke et al. ( 1998 1. 

the j -th surface element. The element is irradiated by each sub¬ 
source with index i under its own angle of incidence //, with 
fj, - cos § j. The irradiation parameter xj for the /-th element is 
given by the first summation term 


with dj the distance between the irradiation source and the /-th 
surface element as depicted in the central graph in Fig. 5. Since 
each surface element is characterized by uniquely defined val¬ 
ues of Xj and pi /, /) can readily be calculated using the spectral 
intensities x j>Pj) from the database. 

The situation differs for an extended source as depicted in 
Fig. 5. The right-hand graph shows a schematic and the left- 
hand graphs the actual model in which the extended source is 
approximated by a set of N subsources of luminosity L\ = L/N 
with N - 100. Because of the curvature of the white dwarf 
surface, only a subset tij < N contributes to the irradiation of 


F in,j _ 1 V 2 ' L, _ (p)j yA Li 

F w d F wd Andr^' F w(i And^ 


( 18 ) 


A complication arises from the fact that the /r, in the sum cover 
a finite range preventing the direct use of the database spectra. 
The exact solution would require us to calculate for each sin¬ 
gle surface element a custom-tailored atmosphere model taking 
into account the irradiation fluxes from all subsources visible 
to that surface element, which would be extremely time con¬ 
suming. For simplicity, we have opted, therefore, to replace the 
actual distribution of //, for the /-th surface element by a flux 
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weighted mean (ji)j, as defined by the right-hand equality in 
Eq. (18), and use the database spectra I,\{6j,Xj, to calcu¬ 
late the observed spectral flux as 

1 m 

h = 757 X A > cos °j Ia( 9 J’ x J’ V>j)- < 19 ) 

j =1 

The radial temperature profiles of the spot created by irradi¬ 
ation from a point source and an extended source differ substan¬ 
tially. As an example, we compare in Fig. |6]the temperature dis¬ 
tributions over the polar cap of a white dwarf irradiated by an 
isotropically emitting point source of luminosity L at a height 
h = 0.1 R vd (solid curves) and by an extended source of the 
same total luminosity represented by an array of N = 100 sub¬ 
sources at the same height and distributed over an area with a 
radius r = 0.3 R w d (dashed curves). We use /? w d = 8.6 x 10 s cm 
appropriate to a 0.6 M 0 white dwarf. The upper and lower set 
of curves are for L = 10 33 erg s 1 and L = 10 32 erg s -1 , respec¬ 
tively. The abscissa is the angular separation a from the center 
of the source, as defined in Fig.|3(center graph). The extended 
source is seen to heat a larger area to a lower peak temper¬ 
ature. These temperature distributions do not account for cy¬ 
clotron beaming, which would broaden them, nor for contribu¬ 
tions from the stream or the X-ray source, which would narrow 
them. 


5. Application to HST/GHRS observations of 
AM Her 

AM Her, the bright prototype of polars, oscillates irregu¬ 
larly between states of high and low (or possibly zero) accre¬ 
tion on time scales of months to years and represents, there¬ 
fore, an ideal laboratory in which to study the irradiation- 
heated white dwarf. Since the discovery that the far ultravi- 
olet continuum of AM Her is dominated by the white dwarf 
jHeise & Verbuntfl (188th numerous spectroscopic studies have 
added to our knowledge of the response of its atmosphere to 
accretion and accretion-induced irradiation. In the low state, 
a deep Lya absorption line is observed at all orbital (= ro¬ 
tational) phases, which bec omes largely filled-up in the high 


state JcLinsicke et all l 1 995l 1 1 9981 20051). The same holds for 


the higher Lyman lines dMauche & Ravmondl ll998l). The lack 


of absorption is illustrated by the conclusion of Greeley et al.l 
ll999t) that the continuum is surprisingly well represented by a 
blackbody. Our calculations add to an understanding of the at¬ 
mospheric structure of the white dwarf in AM Her and provide 
an internally consistent explanation for the lack of deep Lyman 
absorption lines. Since AM Her is typical of the class of polars, 
they add to our understanding of polars, in general. 

For our specific model calculations, we use a distance of 
AM Her D = 79 pc lThorstensenll2003l) . a radius of the white 
dwarf R w[ t = 8.6 x_ 10 8 cm appropriate for a mass of 0.6 M 0 
IGansicke et alj 119951) . and an inclination i = 40°, which is 
within the range of i = 35° - 60° reported for AM Her 
IGansicke et al.ll200ll) . We add interstellar absorption in the 
cores of the Lyman lines using an interstellar column density 
of atomic hydrogen, Ah = (3 ± 1.5) x 10 19 H-atomscm -2 , 
which reproduces the width of the narrow absorption core ob¬ 
served in the Lya line of AM Her ( iGansicke et akifTohii) - The 


continuum spectral shape is not affected noticeably by the ob¬ 
servational uncertainty in the column density. The height and 
radius of the irradiation source are taken as h = 0.1 /? W( j and 
r = 0.3 A\v(i, respectively, which yield the observed spot size of 
~ 10% of the white dwarf surface area (Gansicke et al. 1995; 
iMauche & Ravmondll 1 998t IGansicke et alJl2005h . With these 
assumptions, we calculated the observed spectrum/) as a func¬ 
tion of orbital phase <p. The two left-hand graphs in Fig. 5 illus¬ 
trate the aspect of the white dwarf at two orbital phases, <p = 0 
(best visibility of the spot) and cf> = 0.5 (spot near the horizon, 
source located above the white dwarf). 

Fig. 7 shows the January 1997 high state Hubble Space 
Telescope GHRS spectra of AM Her _at orbital phases 
/ = 0.1 and cf> = 0.5 from IGansicke et al J (Il998l) as 
grey curves (compare the illustrations in Fig. 5). The spec¬ 
tra display strong Doppler-broadened asymmetric emission 
lines of Ciii/11176, Simdl207, LyaT1216, NvT1239,1243, 
Sini/11299,1303, Cn/Cif/11335, and SiivT1394,1403, which 
originate in the accretion stream. We estimate the contin- 
uum flux asso ciated with the stream emission as follows. 
IGansicke el al.j 120051) have compared the FUSE spectra of 
AM Her at different accretion levels and find a wavelength- 
independent stream contribution to the continuum flux for the 
June 2000 high state of 0.4 x 10~ 13 erg cm -2 s -1 A -1 . Converted 
to the brighter January 1997 high state, this contribution is 
about 0.6 x 10~ 13 ergcm _2 s _1 A _I . Pried horskv et. al.l (ll 978l 
their Fig. 1) observed this component with U-band and V-band 
fluxes of 0.6 and 0.2 x 10~ 13 ergcirT 2 s -1 A -1 , respectively, and 
noted that it decreases in the infrared as l/A 2 . Fig. 7 shows 
the calculated spectra at / = 0.1 and d> = 0.5 with the stream 
contribution (dashed line) added. The resulting model spec¬ 
tra are seen to provide an adequate fit to the observed con¬ 
tinuum specta at both orbital phases simultaneously (Fig. 7, 
solid curves). The eye fit fixes the fractional luminosity inter¬ 
cepted and reemitted by the extended spot to 2.3 x 10 32 erg s ', 
of the order of 10% of the bolometric luminosity. At a height 
of 0.1 A\ v( i, the white dwarf subtends a solid angle of 1 ,2n and 
the implied total luminosity of the irradiation source would be 
7.7 x 10 32 erg s ', but this number may be misleading because 
it depends on the actual height of the source and its angular 
emission characteristics. Presently, we prefer simply to com¬ 
pare the integrated model flux for <f> = 0.1 at Earth with the 
observed fluxes of the potential irradiation sources, which con¬ 
tribute to the heating of the polar cap. The model flux longward 
of 912 A at (f> — 0 is 3.3 x 10” 10 erg cm _2 s' 1 , but an additional 
6.0 x 10 _1 °ergcm _2 s _1 or 65% of the total is hidden in the 
Lyman continuum. This large fraction of the Lyman contin¬ 
uum flux has not so far been considered in the energy balace 
of AM Her, but its existence is supported by the observed pro¬ 
nounced Balmer and Paschen jumps in emission in AM Her’s 
high state (iSchachter IHDtmll). The fluxes of the irradiation 
sources at Earth appropriate for the January 1997 high s_tate in¬ 
clude a cyclotron flux of 2.0 x 10 _10 ergcm _2 s _l iBailev et alJ 
I l'OHSl see also Gansicke et al. 1995), a stream contribution for 
A > 912 A of 3.3 x 10 _1 °ergcm _2 s _1 plus probably a substan¬ 
tial Lyman continuum flux, and a hard X-ray flux converted 
from the April 1991 high state value of 1.3 x 10 10 erg cnfV 1 
iGansicke et all 19951) to 1.6xlO _10 ergcm _2 s _I , assuming that 
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the conversion factor is the same as in the FUV. The cyclotron 
flux complemented by fractions of the stream and hard X-ray 
components can account for the observed flux of the heated 
polar cap longward of 912 A. The predicted Lyman contin¬ 
uum flux from the polar cap requires an additional source. 
Candidates are the Lyman continuum emission from the lower 
section of the optically translucent accretion stream and the 
potential soft X-ray emission from a spray of accreted matter 
lifted above the photosphere of the white dwarf as a reaction 
to the drastic fluctuations in ram pressure which occur in non¬ 
stationary accretion. It is clear that the hard X-ray source can 
not play a major role because of its limited luminosity and the 
low conversion efficiency of hard X-rays, as addressed in the 
next section. 


6. Discussion 


We have quantitatively reproduced the 1150 - 1450 A contin¬ 
uum spectra of AM Her in its high state of accretion by a model 
that considers a large polar cap heated by irradiation from the 
accretion region. Irradiation strongly affects the temperature 
stratification of the white dwarf atmosphere in the heated polar 
cap. Strong heating leads to an inversion of the radial temper¬ 
ature profile in the atmosphere, causes the Lyman absorption 
lines to fill up completely, and produces substantial emission in 
the Lyman continuum. The known irradiation sources can ac¬ 
count for the observed reprocessed flux longward of the Lyman 
limit, but an additonal source is needed to account for the pre¬ 
dicted Lyman continuum flux. The implied luminosity of the 
polar cap is larger than considered previously dClansicke et al. 
1995 1 and stresses the point that the reprocessed flux emerges 
primarily in the FUV and not in the soft X-ray regim e as pre¬ 
dicted by the early simple m odels of polars jKinti & L asota 

ll979HLamb & Masters! 1979|). _ _ 

Since X-ray heating ( van Teeseling et mum has not 


been considered here it is appropriate to discuss the poten¬ 
tial contribution from that source in some detail. The hori¬ 
zontal temperature profile of the spot vs. distance from the 
spot center and the vertical atmospheric temperature profile 
at a given position may both be expected to differ for ir¬ 
radiation with either optical/infrared radiation or hard X-ray 
bremsstrahlung. Cyclotron radiation will produce a larger spot 
than bremsstrahlung because it originates higher up in the 
post-shock region and its beaming properties allow it to reach 
regions further away from the spot center. This difference 
is strenghtened by the lower shock heights of the denser, 
bremsstrahlung-dominated shocks (Eq. 14) and by the fact 
that bremsstrahlung is emitted isotropically and irradiation is 
thereby concentrated closer to the spot center. For a given sur¬ 
face element, one can judge the efficiency of heating by cy¬ 
clotron radiation and by bremsstrahlung as follows. As seen 
from Fig. 3 (right panel), cyclotron heating is most intense at 
Rosseland optical depths r r0S s ~ 0.01 but reaches down to the 
photosphere and beyond. Hence, cyclotron radiation is com¬ 
pletely reprocessed fairly high in the atmosphere. Estimating 
the efficiency of bremsstrahlung heating requires us to consider 
the competition between photoabsorption and Compton scat¬ 
tering. Hydrogen is 99% ionized at the relevant levels of our 


atmospheres, which still allows atomic transitions to determine 
T r0S s and causes the electron scattering optical depth to become 
r es ^ Tress- F° r the typical lOkeV bremsstrahlung photons, on 
the other hand, photoabsorption is negligible with r p h <sc r es . 
Hence, bremsstrahlung photons will penetrate far below the 
photosphere and are preferentially (back)scattered rather than 
photoabsorbed with the consequence of ineffective atmospheric 
heating Ivan Teeseling et al.ll994l . Heating by bremsstrahlung 
may be important, however, if a finite metalicity of the atmo¬ 
sphere drastically increases the absorption cross section. Since 
the metalicity of the atmosphere will be enhanced in the in¬ 
ner spot where accretion occurs. X-ray heating may dominate 
there, while cyclotron heating is expected to dominate further 
away from the spot center. While a more detailed model should 
take both irradiation sources into account, we presently lack in¬ 
formation on parameters like the mass inflow rate per unit area 
and the ensuing metalicity as functions of distance from the 
spot center. 

The heating of the extended polar cap from the post-shock 
source tends to zero for vanishing shock height because then all 
downward directed emission is intercepted by the atmosphere 
within the proper accretion spot. Disregarding irradiation by 
emission from the pre-shock stream for the moment, our model 
predicts that high-field polars, in which all cyclotron and X-ray 
emission occurs in low-lying shocks (see Eq. 14), should lack 
strongly irradiation-heated polar caps. This is, in fact, the case 
for the 200 MG system AR UMa in its low state with some rem¬ 
nant cyclotron emission and an amp litude of the orbital mod¬ 
ulation at 1400A of only 5% iGansicke et alJEoPll) . An inves¬ 
tigation of the temperature distributions of a larger sample of 
white dwarfs in polars would shed light on some poorly under¬ 
stood aspects of accretion physics. Such a study would have to 
consider cooling of the deeply heated polar cap as well as the 
relevent heating processes. Note that accreting white dwarfs 
are generally hotter than single white dwarfs of the same age 
dArauio-B etancor et. al J20 05l)as a result of compressional heat- 
ing JTownslev&Biidstenl20031) . but this process is rooted in 
the entire non-degenerate envelope (and to some extent in the 
core) and heats the entire star and not just the polar cap where 
accretion takes place. 

Finally, a source of uncertainty in the calculated spec¬ 
tra warrants mentioning. Our atmospheric model disregards 
the compression of the outer atmosphere by accretion and 
is applicable strictly only outside the region of the extended 
source, assumed to be fed by accretion. The bombardment of 
the atmosphere by charged particles adds an external pres¬ 
sure and affects the structure of the upper atmosphere where 
the emission lines originate. The physics of bombarded atmo¬ 
spheres is_complicated and has been treated to a certain ex¬ 
tent bv IWoelk & Beuermanni 11 9921 ). A simple estimate sug¬ 
gests that the effect may not be negligible. For a mass flow 
rate m = 10 4 g cm 2 s l , which is below present detection 
techniques, the ram pressure P mm = mvg, with vg the free- 
fall velocity, amounts to a few percent of the photospheric 
pressure and severely compresses the atmosphere only outside 
Uoss — 3 X 10 4 . A mass flow rate ten times higher, however, 
reaches down to r ross - 0.05. Under these circumstances, a 
more detailed study of the line profiles and the strength of the 
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Lyman continuum in irradiated atmospheres of mCVs, includ¬ 
ing NLTE effects, requires a substantially larger effort. 
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